The Journal of Pathology

Primary sclerosing cholangitis (PSC) is a rare, progressive liver disease characterized by biliary inflammation and bile duct strictures and no approved medical therapy. Despite its clinical severity, the pathological mechanisms underlying PSC remain poorly understood, largely due to early diagnostic challenges. Here we provide complementary evidence in human PSC samples, transcriptomic data, mouse models, and 3D cholangiocyte cultures that underscore the importance of hyaluronic acid (HA) and its cognate receptor CD44 in PSC pathogenesis. HA is a glycosaminoglycan abundant in the extracellular matrix in inflammatory disorders, yet its role in PSC has not been well characterized. We demonstrate that in early-stage PSC, cholangiocytes aberrantly produce high molecular weight HA that accumulates in the peribiliary matrix, increasing local tissue stiffness. This mechanical signal is transduced by a CD44/Integrin {beta}1 receptor complex in cholangiocytes, driving cell proliferation, YAP mechanosignaling, pro-inflammatory cytokine production with a transition to a ductular reactive phenotype. CD44 knockdown in cholangiocyte cell lines and mouse models significantly lowered stiffness, and attenuated inflammation. Together, these findings reveal a mechano-inflammatory axis in which HA-driven matrix stiffening perpetuates biliary inflammation and disease progression, identifying HA targeting and CD44 as promising therapeutic strategies. One Sentence SummaryHyaluronan and CD44 mediate matrix changes and progressive fibrosis in primary sclerosing cholangitis.

Background: Chronic inflammation is the most common histological feature in Benign Prostatic Hyperplasia (BPH), and T cells are a key component of immune infiltrate. Advanced BPH is commonly associated with the formation of nodules, but it remains unclear whether a link exists among T cell infiltration, nodular development, and BPH progression. Using a Toxoplasma gondii (T. gondii) model and human specimens, we characterize the subtypes of T cells present during prostatic hyperplasia and their association with nodular development of the prostate. Methods: Male CBA/j mice were intraperitoneally infected with T. gondii parasites, and flow cytometry was performed on the prostate to quantify the number of CD4+ and CD8+ T cells. Histology was used to score microglandular hyperplasia (MGH), and immunofluorescence was used to quantify and examine the locality of CD4+ and CD8+ T cells and compared that to human BPH tissue. Results: We found that infecting male mice with T. gondii resulted in an increase of both CD4+ and CD8+ T cells in the prostate acutely and that CD8+ cells remained sustained at chronically. We also established the presence of glandular nodule formation at this timepoint through hematoxylin and eosin (H&E) staining. Immunofluorescence revealed that CD8+ cells were found proximal to forming glandular nodules relative to non-nodular glands. We also found more CD8+ cells localized to non-nodular glands in nodular BPH tissue versus non-nodular BPH tissue. Finally, we discovered a higher prevalence of CD8+ cells in T. gondii IgG+ patients than in IgG- patients. All T. gondii IgG+ patients exhibited nodular BPH, whereas all but one IgG- patient exhibited non-nodular BPH. Conclusions: This study is the first to investigate the presence and location of CD4+ and CD8+ T cells within nodular and non-nodular BPH glands. We found an association of the presence of CD8+ T cells with nodular progression. This association held true in human prostate tissue. Translationally, CD8+ T cells may enhance nodular BPH progression, and T. gondii infection may promote this CD8+ T cell-mediated response.

Epidermodysplasia verruciformis (EV) is characterized by susceptibility to cutaneous {beta}-human papillomaviruses, persistent viral lesions and an increased risk of cutaneous squamous cell carcinoma (cSCC). We investigated whether expression of Toll-like receptors 3, 4, 5 and 9 (TLR3, TLR4, TLR5 and TLR9) differs across histological states within EV. Immunoreactivity was quantified by digital H-score using QuPath in 97 specimens from 28 patients, comprising 13 lesion-free skin samples, 15 flat warts and 69 cSCCs. Prespecified comparisons evaluated flat wart versus lesion-free skin, cSCC versus non-neoplastic tissue and invasive versus in situ cSCC, with invasive tumour differentiation assessed secondarily. Patient-clustered generalized estimating equations and false-discovery-rate (FDR) correction were used for inference. After correction for multiple testing, TLR9 H-scores were higher in flat warts than in lesion-free skin ({beta}=60.58; q=0.0332). Compared with non-neoplastic tissue, cSCC showed lower TLR3 ({beta}=-81.40; q=0.00745), higher TLR4 ({beta}=55.59; q=0.04793) and lower TLR9 expression ({beta}=-34.80; q=0.03412). TLR5 showed a nominal increase in cSCC that did not meet the FDR threshold. None of the four receptors differed significantly between invasive and in situ cSCC, and no differentiation trend remained significant after FDR correction. These findings show receptor-specific differences in TLR expression across histological states of EV, particularly involving TLR3, TLR4 and TLR9, although the diagnostic and prognostic significance of these tissue patterns remains to be established.

Intraductal papillary mucinous neoplasms (IPMNs) are cystic precursors of pancreatic ductal adenocarcinoma undergoing dynamic epithelial and stromal remodeling during progression. In this study, we integrated computational pathology, cyclic immunofluorescence, and Xenium spatial transcriptomics to define the spatial organization of the IPMN microenvironment. We queried cell-based histopathological features and identified stromal patterns strongly associated with epithelium proximity. Cyclic immunofluorescence and Xenium revealed a subepithelial gradient extending from myofibroblasts to inflammatory fibroblasts. Xenium further identified a distinct subepithelial NOTCH1-NESTIN-expressing fibroblast subset in myofibroblasts that expands with tumor progression. By recapitulating the spatial remodeling event using human pancreatic tumor cells and primary myofibroblast co-culture models, we demonstrated the capacity of pancreatic tumor cells to induce a NESTIN-expressing state in neighboring myofibroblasts. These findings reveal new dynamic epithelial-stromal interactions through subepithelial NESTIN+ fibroblasts associated with pancreatic tumorigenesis and immediate microenvironment remodeling by active tumor cells, highlighting the power of computational pathology integrated with spatial transcriptomics.

Muscle-invasive bladder cancer (MIBC) carries high risks of relapse and death despite contemporary care. In this study we aimed to evaluate whether the pretreatment immune microenvironment varied with HER2 status and if the differences influenced the response to neoadjuvant chemotherapy (NAC) and survival. We conducted a retrospective study of 62 treatment-naiive MIBC patients diagnosed on pretreatment transurethral resection of bladder tumor (TURBT) specimens collected between 2008 and 2024. All patients received NAC followed by radical cystectomy. HER2 immunohistochemistry (IHC) was scored 0/1+/2+/3+ using gastroesophageal criteria. PD-L1 was assessed by combined positive score (CPS). CD3, CD20, and CD68 cell densities were quantified by digital image analysis. Overall survival was measured from the date of TURBT. HER2 IHC scores were 2+/3+ in 28 of 62 tumors (45%) and 0/1+ in 34 (55%). All seven IHC 2+ tumors lacked ERBB2 amplification by fluorescence in situ hybridization (FISH). NAC response occurred in 11 of 34 tumors with HER2 IHC 0/1+ (32%) and 9 of 28 tumors with HER2 IHC 2+/3+ (32%) (P = 1.00). PD-L1 expression was not associated with HER2 status. Immune-cell analyses showed trends toward higher CD3+ T-cell density, CD20+ B-cell density, and TM ratio in NAC responders, but no significant differences by HER2 status or combined PD-L1/HER2 subgroup status. After Bonferroni correction for multiple comparisons, none of the immune-cell metrics remained statistically significant. Overall survival did not differ between the HER2 IHC groups (P = 0.71). These findings indicate that HER2 IHC 2+/3+ expression is present in a substantial subset of MIBC but is not associated with NAC response or overall survival in this cohort. Further study is warranted to define the clinical relevance of MIBC with HER2 IHC 2+/3+ and its potential role in HER2-directed therapeutic strategies.

Under homeostatic conditions, the cornea is avascular and contains few immune cells, but this changes rapidly following injury. Although the long-term consequences of corneal damage are well characterized, the earliest vascular and immune responses remain poorly understood. Here, we used a murine corneal alkali-burn model to examine limbal vascular activation and leukocyte recruitment immediately and at 2, 6, and 24 hours after injury. Limbal blood vessels underwent immediate dilation; however, vascular leakage into the corneal stroma occurred only in males. Lymphatic capillaries rapidly formed directed extensions toward the injury without significantly increasing their total vascular area, with males exhibiting longer extensions than females. Fluorescent dextran uptake provided evidence that these lymphatic vessels were functionally engaged in early tracer drainage. Despite pronounced vascular activation, early recruitment of neutrophils, monocytes, dendritic cells, macrophages, T cells, B cells, and natural killer cells remained limited. Thus, limbal blood and lymphatic vessels initiate the earliest response to corneal alkali injury before substantial leukocyte infiltration. These findings reveal sex-dependent differences in vascular permeability and lymphatic remodeling and identify the limbal vasculature as an early regulator of corneal inflammation and tissue repair.

Abstract Background: Intraductal papillary mucinous neoplasms (IPMNs) are recognized as precursor lesions to pancreatic ductal adenocarcinoma (PDAC). However, the molecular programs underlying progression from low-grade dysplasia to advanced disease remain incompletely characterized. Herein, we performed an integrated plasma and tissue-proteomic analyses coupled with spatial and single-cell transcriptomics to identify biologically coherent remodeling programs reflected in circulation that distinguish IPMN by dysplasia grade and invasive disease. Methods: Using the O-link proximity extension assay platform, a panel of 1,104 proteins were quantified in plasma samples collected from patients with low-grade (LG) IPMN (n=30), high-grade (HG) IPMN with or without associated PDAC (IPMN/PDAC; n=40) and PDAC without IPMN (n=8). Predictive performance of individual biomarkers were assessed; likelihood ratio testing was performed to identify protein biomarkers that were complementarity with CA19-9 for risk of malignancy of IPMN. Findings were intersected with available spatial (N= 13) and single-cell (N= 6) transcriptomic datasets of IPMN tissues as well as mass spectrometry-based proteomic profiles of an independent set of resected human IPMN tissues (N= 9). Results: A total of 28, 43, and 35 circulating proteins were found to be differential in HG, IPMN/PDAC, and HG + IPMN/PDAC cases compared to LG IPMN. Among differential proteins were known PDAC-associated markers CEACAM5, CTRC, and REG3A as well as several biomarkers reflecting cytoskeletal and extracellular matrix remodeling and inflammatory processes. Focusing on cytoskeletal and ECM-related proteins and using likelihood ratio testing, an OR rule considering CA19-9, BGN, and ITGB1BP1 achieved overall sensitivity of 48.7% for HG + IPMN/PDAC, including 38.1% sensitivity for HG IPMN, at an overall specificity of 90%, which was improved compared to that of CA19-9 alone (overall sensitivity of 28.2%; McNemar Exact test 1-sided p-value: 0.011). Integrated proteomic and spatial transcriptomic datasets of IPMN tissues revealed coordinated alterations cytoskeletal and ECM remodeling and elevated matrix stiffness as prominent features associated with IPMN/PDAC, which paralleled concordant increases in BGN and ITGB1BP1. Cell-type of origin analyses based on spatial and single-cell data further revealed fibroblasts and myeloid cells as primary contributors to expression levels of BGN whereas ITGB1BP1 was primarily expressed in neoplastic epithelium. Conclusion: Advanced IPMN dysplasia and invasive disease are characterized by coordinated tissue remodeling programs that are systemically reflected in circulating proteomic profiles. Blood-based biomarkers identified through our study, such as BGN and ITB1BP1, have potential to improve upon CA19-9 for risk stratification of IPMN to better guide clinical management.

Objective: The objective of this study was to detect genetic factors associated with dermatochalasis using a genome-wide association study (GWAS) across three large cohorts. Design: GWAS meta-analysis Participants: A total of 13,200 dermatochalasis cases and 962,513 controls were included. Methods: A GWAS meta-analysis of dermatochalasis combining data from the FinnGen, the Estonian Biobank and the UK Biobank was conducted. We also performed colocalization analyses, a phenome-wide association study and age-at-onset analysis, and assessed genetic correlations with various diseases and traits. Main outcome measures: Identification of genetic variants associated with dermatochalasis. Results: We identified 18 loci associated with dermatochalasis at genome-wide significance, 16 of which were novel. Most of these loci had genes involved in skin biology and cutaneous diseases, such as the genes encoding elastin (ELN) and Latent TGF-{beta} binding protein 1 (LTBP1). Phenome-wide association study revealed previous associations with morphology-related traits, while genetic correlation analysis highlighted multiple genetic correlations, especially with smoking and pain. Conclusions: We detected 18 genetic loci associated with dermatochalasis, characterized these loci in detail and demonstrated their relevance in skin biology and related processes. These findings give novel information on the genetic background of dermatochalasis and provide a solid basis for further research.

Genetic variants in PROM1 are associated with inherited blindness, but clinical phenotypes vary widely with currently no consensus on disease expectations or predicted patient outcomes. To address this issue, we performed chi-square correlation analysis on 190 published pathogenic and likely pathogenic variants from ClinVar Miner and the Human Gene Mutation Database. We identified a statistically significant increase in observed cases of dominantly inherited missense variants, which were associated with milder forms of macular dystrophy. Recessive inheritance of truncating frameshift variants was associated with more severe forms of inherited retinal degeneration such as rod-cone dystrophy. We also present case studies from 7 patients with phenotypes that are unique, unreported, or do not match the reported phenotype in online databases. We report two patients with VUS that have clinical phenotypes and dominantly inherited missense variants, one patient with a long natural history at our clinic, two patients with the same variant, but a different clinical diagnosis, and two patients with retinal deposits comprising subretinal drusenoid deposits (SDD) and soft drusen.

Background: Pre-operative triage of adnexal masses with serum CA-125, HE4, and ultrasound (O-RADS) has limited accuracy, contributing to unnecessary surgery. We develop and validate a first of its kind, biomarker panel-free plasma proteomic classifier that distinguishes malignant from benign adnexal masses from a single blood draw, operating directly on discovery-mode mass spectrometry. Methods: In this multicenter prospective observational study, plasma from women with adnexal masses and from healthy controls at a large urban academic center (discovery and four validation cohorts) and an external academic center from another Urban-suburban area (external validation) was analyzed by data-independent-acquisition mass spectrometry (Orbitrap Astral), quantifying >1,000 proteins per sample. Cancerous (malignant, borderline adnexal mass, or metastasis to ovary) versus benign adnexal masses were classified with ADAPT-MS (Adaptive Diagnostic Architecture for Personalized Testing by Mass Spectrometry), which retrains an ElasticNet model for each sample based on the proteins measured. The primary outcome was discrimination (area under the receiver operating characteristic curve [AUC]). Pre-operative O-RADS scores and serum CA-125 were compared and combined with proteomics on matched subsets. Findings: Using a discovery cohort (n=1062), an ADAPT-MS classifier, validated across five independent cohorts (n=623), achieved a summary AUC of 0.835 (cohort range 0.779-0.904). Among patients with pre-operative O-RADS (n=231), proteomics (AUC 0.852) was comparable to O-RADS alone (0.863) and combining them raised discrimination to 0.922. Proteomics outperformed serum CA-125 (0.915 versus 0.777; n=125) and proteome-derived ROMA-like and two-marker surrogates. Biologically, the classification-relevant proteins are enriched for host response proteins rather than tumor cell associated protein changes, comprising ECM remodeling, complement and innate immune response pathways. This explains good detection of tumor with heterogenous pathology and also metastasis to ovary. Interpretation: A single plasma proteomic measurement interpreted using adaptive machine learning discriminates malignant from benign adnexal masses with good internal and external validation, and is complementary to ultrasound, warranting prospective evaluation as a panel-free diagnostic adjunct to imaging.

Breast cancer (BC) shows greatest sexual diversity. Increased diagnosis and poorer outcomes in men highlights the need to better understand its biology. We hypothesised that cancer-associated fibroblasts (CAFs), the most abundant cell type in the tumour microenvironment, might define sex-related differences. Using phenotypically matched male and female CAFs generated from breast cancer tissues, we demonstrate distinct transcriptional programmes and functional behaviours associated with extracellular matrix remodelling, cell adhesion, migration and vascular development. Compared to CAFs generated from females BC, those from males generated denser, more complex matrices promoting stronger tumour and endothelial cell adhesion, vascular growth, but less organised capillary network formation. Findings reveal fundamental sex-related variations in CAF phenotype and biology in BC. These findings highlight the need to integrate biological sex into precision oncology to identify opportunities for sex-specific therapeutic strategies in BC.

Background & Aims: Cholangiocarcinoma (CCA) is an aggressive malignancy with poor five-year survival. Improving outcomes requires preclinical models that faithfully recapitulate the native tumour microenvironment. Precision-Cut Tissue Slices (PCTS) retain in-vivo architecture and represent a promising candidate. We aimed to establish a patient-derived PCTS model of CCA, characterise its response to ex-vivo culture, and assess its use as a platform for therapeutic testing. Methods: PCTS were generated from 25 patients undergoing curative-intent resection for presumed CCA (2022-2024) and cultured for up to 15 days. Viability and histological architecture were assessed serially. Quantitative proteomics (SWATH-DIA) was performed at Days 0, 3, 7 and 15. A subset were treated with staurosporine or clinically relevant chemotherapy (5'-deoxy-5-fluorouridine, gemcitabine {+/-} cisplatin). Results: PCTS viability was maintained to Day 15, with no significant reduction from baseline. PCTS retain tumour architecture, cytokeratin-19-positivity and resident CD3+/CD68+ immune cells during culture. Proteomic profiling quantified 4,578 proteins and identified a staged response, with differentially abundant proteins increasing from 8 (Day 3) to 127 (Day 7) and 180 (Day 15), stabilising thereafter. This comprised of an early loss of inflammatory and stromal proteins, sustained cellular stress responses and metabolic reprogramming. CCA subtype-specific proteomic signatures were preserved throughout culture. Staurosporine produced dose-dependent cytotoxicity whilst chemotherapy responses were variable. Conclusions: Patient-derived CCA PCTS maintain viability, tumour architecture and subtype identity for 15 days while undergoing a defined proteomic response to culture, and retain pharmacological responsiveness,

Individuals with fibrodysplasia ossificans progressiva (FOP), a rare genetic disorder caused by mutations in the bone morphogenetic protein receptor ACVR1 (also known as ALK2), experience progressive and severely debilitating endochondral heterotopic ossification (HO). Previous studies have identified fibro-adipogenic progenitors (FAPs) as a major cellular source of HO in FOP mouse models. Here, we used mRNA sequencing of FAPs isolated from injured muscle of FOP mice to capture the earliest gene expression changes associated with the transformation of muscle tissues into bone. Whereas the transcriptomes of mutant and wild-type FAPs were indistinguishable prior to injury and at 1 day post injury, dramatic dysregulation of gene expression was observed at 3 and 5 days after injury. Surprisingly, gene expression changes suggested early specification of osteochondral lineages by 3 days post injury, approximately 3 to 4 days before the appearance of histologically identifiable cartilage. We identified dysregulation of genes encoding components of major signaling pathways, cell cycle regulators, extracellular matrix components, and hypoxia response mediators, and showed that these gene expression changes are temporally associated with abnormal population growth of mutant FAPs. We further demonstrated that treatment of FOP mice with either an anti-activin A monoclonal antibody or the ACVR1 kinase inhibitor BLU-782 (IPN60130; fidrisertib) normalized the expression of the majority of dysregulated FAP genes. However, the relative efficacy of the drugs differed, and distinct sets of overlapping genes were normalized by the treatments. The substantial transcriptional changes caused by these drugs in FAPs from normally regenerating muscle were also enumerated and described. These findings define early transcriptional events associated with pathogenic FAP reprogramming and provide a molecular framework for understanding the shared and distinct effects of activin A and ACVR1 inhibition.

Background: Colorectal cancer liver metastasis is accompanied by coordinated stromal remodelling, but which fibroblast and malignant epithelial states are coupled, and whether any such coupling holds outside the dataset it was found in, is rarely tested. Methods: Paired primary tumours and liver metastases from six patients were analysed by single-cell RNA sequencing with copy-number-based identification of malignant cells, trajectory inference, cell-cell communication and ligand-target modelling, and spatial transcriptomics of three metastatic sections. Protein staining was taken from an independent tissue resource. Every step was then repeated in a second, independent paired cohort of five patients and two further metastatic sections, with cell identities fixed by transfer from the discovery data rather than refitted. Results: An RGS5-positive fibroblast/pericyte-like state (F2) and a WNT-like malignant epithelial state (C2) were both enriched in metastases and occupied late trajectory positions. Communication analysis nominated a directional F2-to-C2 programme with transforming growth factor beta as the highest-ranked pathway and its receptor TGFBR2 enriched in C2. Spatial data placed the two states in the same neighbourhoods. Hub genes of the programme were TPT1, SLC1A5, SOX4 and TSC22D1, and protein staining was higher in tumour than normal colon for the three with available data. In the independent cohort the same states reappeared: F2 was enriched in metastases (Ro/E 1.22 versus 0.52; 32% versus 14% of fibroblasts), C2-like cells were the most metastasis-enriched epithelial state (1.10 versus 0.86), 92% of the cells that cohort's own authors called MCAM-positive fibroblasts mapped to F2, and F2 was again the strongest TGFB sender to C2-like epithelium at both sites. Conclusions: A directional TGFB axis from an RGS5-positive fibroblast state to a WNT-like epithelial state is a reproducible feature of colorectal liver metastasis rather than a property of one dataset. Its spatial expression varies between sections, which constrains how tightly the two states can be said to be co-localised.

Background and Objectives: Facioscapulohumeral muscular dystrophy (FSHD) is characterised by asynchronous muscle degeneration and marked pathological heterogeneity. While conventional MRI can detect fatty replacement and oedema-like abnormalities, it provides limited insight into other relevant biological processes such as extracellular matrix remodelling and tissue perfusion. This study investigated whether multispectral optoacoustic tomography (MSOT), an emerging imaging technique that non-invasively interrogates tissue composition through endogenous optical absorbers, can provide complementary molecular information on muscle involvement in FSHD. Methods: MSOT was performed on 568 muscles from 15 genetically confirmed FSHD patients and 9 healthy volunteers using a standardized protocol. The analysis included single-wavelength signal intensity and spectrally unmixed signals related to lipid, collagen, deoxygenated haemoglobin, oxygenated haemoglobin, total haemoglobin, and tissue oxygen saturation. Patients also underwent conventional muscle MRI together with clinical strength and disease severity assessments. Feasibility, reproducibility, discrimination between patients and controls, and associations with MRI and clinical measures were evaluated. Results: MSOT acquisition was well tolerated and highly reproducible (intra- and inter-rater ICC = 0.88). FSHD muscles displayed a distinct optoacoustic profile characterised by reduced signal intensity at wavelengths below 900 nm, increased lipid signal and reduced oxygenated and total haemoglobin signals. Lipid signal discriminated FSHD from control muscles (AUC 0.830, p = 0.008) and increased with MRI-defined fatty replacement (p < 0.0001). Short tau inversion recovery (STIR)-positive muscles on MRI exhibited a distinct molecular phenotype, with higher collagen-associated and lipid signals together with increased deoxygenated and total haemoglobin signals, consistent with active tissue remodelling. Although FSHD muscles classified as normal by conventional imaging showed group-level optoacoustic profiles comparable to controls, 13 of 141 (9%) demonstrated at least one abnormal MSOT-derived parameter. Reduced haemoglobin-related signals correlated with greater global disease severity, whereas increased lipid signal was associated with lower muscle strength. Discussion: MSOT is a feasible and reproducible imaging technique in FSHD that provides complementary molecular information beyond conventional structural imaging. Its ability to detect signatures related to fatty replacement, STIR positivity, functional impairment, and isolated abnormalities in radiologically normal muscles supports further longitudinal studies of MSOT-derived readouts as potential biomarkers of disease activity, progression, and therapeutic response.

When breast cancer invasion begins and how tumor cells breach the basal barrier remain poorly defined. We profiled normal mammary ducts, ductal hyperplasia (DH), ductal carcinoma in situ (DCIS) and invasive ductal carcinoma using spatial transcriptomics, spatial proteomics and five bulk-omics layers, alongside an independent longitudinal lesion cohort. Cross-sectionally, basal/myoepithelial continuity declined most between DH and DCIS, accompanied by extracellular-matrix remodeling and altered fibroblast- and macrophage-associated signaling. EGFR-positive luminal progenitor-like cells were enriched at manually annotated basal discontinuities and were molecularly distinct, nominating a candidate leader-like population without establishing causality. In the longitudinal cohort, expression of GABRG3, TAGLN, MLPH and AZGP1 in initially benign lesions was associated with subsequent ipsilateral malignancy. These findings support a model in which progression-relevant breast tissue remodeling may begin at the DH stage and nominate cellular states and candidate biomarkers for prospective validation in breast cancer risk stratification among patients with DH.

Leiomyosarcoma is an aggressive soft tissue sarcoma that metastasizes preferentially to lung and liver. Because matched primary and metastatic specimens are rarely available, the cellular changes that accompany metastatic progression remain poorly defined. Here we combined single-nucleus RNA sequencing of archival FFPE tissue from four matched primary-metastasis pairs (two lung, two liver metastases; 48,941 nuclei) with imaging-based spatial transcriptomics of an extended tissue microarray cohort of 21 patients (69 cores, 313,218 cells). Malignant cells, not assignable from transcriptomic references alone, were defined by inferred copy number profiles that were concordant with matched DNA methylation arrays and confirmed the clonal origin of each metastasis. Malignant metaprograms were organ-associated, with a stress- and hypoxia-associated program enriched in lung metastases and a proliferative, matrix-remodeling program enriched in liver metastases. Antigen-presenting cancer-associated fibroblasts (apCAFs) were enriched in metastases relative to matched primaries. This shift was reproduced in the spatial cohort (p = 0.001), where apCAFs formed a discrete band at the tumor-organ interface between malignant cells and adjacent parenchyma, whereas they were sparse and randomly distributed in primary tumors. Spatial transcriptomics further resolved the malignant compartment into three histomolecular phenotypes: (i) smooth-muscle-like with myxoid features, (ii) GPNMB with anaplastic features, and (iii) PDGFRB with vessel-rich stroma, each with a concordant immunohistochemical (-SMA, GPNMB, PDGFRB) and morphological correlate. Together, our data provide a matched primary-metastasis reference for leiomyosarcoma, identify apCAFs as a spatially organized, metastasis-associated stromal population, and delineate histomolecular phenotypes recognizable by routine immunohistochemistry, including a GPNMB phenotype potentially amenable to GPNMB-directed therapies.

Background: Pancreatic ductal adenocarcinoma (PDAC) is projected to become the second leading cause of cancer-related mortality, with an increasing incidence and dismal 10% survival at 5 years despite best available care. PDAC is characterized by low tumour cellularity, recurrent driver mutation alterations, and substantial inter and intra-tumoural heterogeneity. The influence of individual genetic features such as sex and ancestry on tumour heterogeneity remains poorly understood. Methods: We assembled a cohort of PDAC across five datasets for a total of 7,527 tumours. Measures of genomic instability were assessed across sex and ancestry on a global level (SNV density and percent genome altered by copy number events), and at a gene level for recurrent events. KRAS variant frequencies, ploidy and mRNA abundance were also evaluated. Ancestry was assigned using ADMIXTURE from germline genotypes where available, and patient reported otherwise. Results: We observed ancestry-specific mutation patterns in global point mutation density, with higher SNV burden in African PDAC tumours and lower SNV burden in Asian PDAC tumours as compared to European tumours. We also identified an ancestry bias in the driver gene, GNAS, in which driver SNVs in GNAS was detected more frequently in those of Asian ancestry, and less frequently in those of African ancestry. No sex- or ancestry-associations were found for KRAS, but individual KRAS variants showed a trend towards ancestry-related variation. Conclusions: These findings demonstrate that ancestry contributes to the genomic heterogeneity of PDAC, shaping both global and gene-specific mutation patterns. This work underscores the importance of including diverse ancestries in PDAC molecular profiling efforts which may inform the development of tailored molecular diagnostics and targeted therapies.

Malignant peripheral nerve sheath tumors (MPNSTs) are aggressive sarcomas. An obstacle to treating MPNSTs is a lack of effective systemic therapies. Although over 70% of human MPNSTs have lost or inactivated the epigenome regulator polycomb repressive complex 2 (PRC2), its activity and contribution to canine PNST progression remain unclear. This study compared canine peripheral nerve sheath tumors (PNSTs) and human MPNSTs across biological and clinical features, including PRC2 activity. Immunohistochemical analysis was performed for a human tissue microarray of 54 neurofibromas and 139 MPNSTs, and 63 canine PNSTs for H3K27me3, a repressive histone mark deposited by intact PRC2, and H3K27ac, which increases globally upon H3K27me3 loss. To understand the genomic alterations present in canine PNSTs, we analyzed tumor mutation burden, copy number alteration, and transcriptomes of eight canine PNST/normal pairs. The results suggested that H3K27me3 loss and associated gain of H3K27ac epigenetically drive human and canine tumors. These findings warrant further studies to evaluate whether these epigenetic deregulations alter similar gene signatures across species.

Background & AimsHaematoxylin and eosin (H&E) staining remains the diagnostic gold standard for solid cancers, including hepatocellular carcinoma, and is increasingly complemented by genomic profiling for precision medicine. Inferring genomic alterations directly from H&E images could streamline testing, but heterogeneity and biases in human training data limit interpretation of genotype-phenotype associations. Here, we aimed to relate histologic to genomic pathology to provide biological explainability for mutation prediction models and assess the impact of germline variation on model performance. MethodsWe analysed 597 murine liver tumours with matched whole-genome sequencing and histopathology (163,835 image tiles; 22.9 million nuclei). Our controlled in vivo design accounted for germline variation, biological sex, and causal mutagen (N-diethylnitrosamine), removing confounding factors present in human cohorts. We trained and evaluated deep learning and supervised machine learning models to predict germline variation and cancer driver alterations from H&E. ResultsModelling accurately predicted germline and somatic alterations from histology, at both locus-specific and genome-wide scales. Quantitative image analysis revealed an unexpected association between Egfr driver mutations and hepatic steatosis, linking genotype to an interpretable morphological phenotype. While model performance declined when applied to tumours from unrepresented genetic backgrounds, this limitation was biologically informative, revealing strain-dependent differences in tumour evolution, notably the prevalence of whole-genome duplication. ConclusionsMachine learning integration of histological and genomic pathology enables accurate, interpretable inference of genetic alterations from H&E, potentially reducing reliance on costly ancillary molecular assays. Our predictions are supported by human-interpretable biological features, addressing concerns around "black-box" technologies. However, caution is required when applying such methods to samples with a genetic background that, even if closely related, is beyond the genetic horizon of training data.